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ABSTRACT  ! 

1 

The  dose  in  rep.  is  calculated  at  interior  points  of  a 
half  space  of  soft  tissue,  for  the  case  of  neutrons  entering  j 

the  tissue  from  a non-capturing  isotropically  scattering  j 

medium.  The  variation  of  mean  free  path  with  energy  is 
neglected,  and,  with  this  approximation,  the  biological  effects 
can  be  shown  to  be  proportional  to  the  energy  flux.  This  is 
then  expressed  as  a function  of  depth  of  tissue  using  the 
Spherical  Harmonics  P7  approximation. 
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1.  Introduction. 


The  object  of  this  report  is  to  calculate  the  biological  effect  on 
the  human  body  when  monochromatic  fast  neutrons  fall  nearly  isotropically 
on  the  surface  of  the  tissue.  The  neutrons  are  supposed  to  travel 
through  a non-capturing  isotropically  scattering  medium  in  direct  contact 
with  the  tissue.  Consequently  neutrons  enter  the  tissue  in  nearly 
isotropic  distribution. 

The  calculation  is  similar  to  that  carried  out  in  T/R  275  where  the 
calculation  of  the  energy  deposition  in  soft  tissue  was  made  for  neutrons 
entering  at  right  angles  to  the  free  surface  (the  albedo  case). 

The  conditions  assumed  in  the  present  report  correspond  to  a person 
leaning  against  the  reflector  of  a pile  and  for  convenience  we  shallrefer 
to  the  non-capturing  medium  as  the  "reflector".  For  most  practical 
considerations,  such  as  neutrons  emerging  from  the  shield  of  a pile, 
neutrons  will  enter  the  body  in  a more  anistropic  manner  than  those 
coming  from  the  reflector  because  the  neutrons  which  have  the  greatest 
chance  of  penetrating  the  shield  are  those  which  have  travelled  the 
least  distance  in  it,  i.e.  those  emerging  from  the  shield  of  the  pile  at 
right  angles  to  it. 

The  calculations  of  the  present  report  and  those  of  T/R  273  cover 
the  two  extreme  cases  of  normal  incidence  and  almost  isotropic  incidence 
and  they  should  be  sufficient  to  assess  the  energy  deposition  in  tissue 
for  the  intermediate  cases. 

The  present  calculations  differ  from  those  reported  in  T/R  273  in 
another  detail.  In  the  present  case  the  body  is  in  contact  with  a 
scattering  medium  so  that  the  neutrons  may  be  scattered  successively 
between  the  body  and  the  reflector.  In  T/R  273  the  body  was  supposed  to 
be  effectively  in  vacuo.  However  the  conditions  of  T/R  273  apply  also 
if  the  body  is  in  direct  contact  with  a non-scattering  absorbing  medium 
and  in  this  sense  the  two  reports  cover  the  two  extreme  cases  of  a pure 
absorber  and  a pure  scatterer. 

Successive  scattering  between  the_body  and  the  reflector  will  of 
course,  be  negligible  if  the  body  is  at  a distance  from  the  reflector 
which  is  large  compared  with  the  dimensions  of  the  body.  The  present 
calculation  is  for  most  purposes  pessimistic  because  it  assumes  direct 
contact. 

The  assumptions  made  in  this  analysis  are  the  same  as  those  in  T/R  273* 

The  energy  transferred  to  protons,  which  is  a measure  of  the 
ionisation  produced  in  a given  volume  element  per  unit  time  is  given  by 


■I  / E 


n(.r,  v)  vdv 
ip(  v) 


where  ip  (v)  is  the  mean  free  path  in  tissue  for  collisicns  with  a proton, 
n(  v)dv  is  the  density  of  neutrons  with  velocity  v in  the  velocity 
interval  dv,  at  the  elementary  volume  whose  position  vector  is_r  and  E is 
the  energy  of  the  neutrons. 

As  in  T/R  273>  • ipt  v)  will  be  assumed  constant.  The  over-estimations 
of  neutron  density  and  mean  free  path  act  in  opposite  directions  in  the 
integral  and,  to  some  extent,  compensate  each  other. 


It  is  true,  also,  that  the  neutrons  which  make  the  greatest  contribution 
to  the  integral  are  those  with  highest  energy.  Hence  the  assumption  of 
constant  mean  free  path  is  justified.  With  this  approximation,  the 
integral  is  proportional  to  the  neutron  energy  flux. 

For  simplification,  the  human  body  will  be  replaced  by  a half  space 
of  soft  tissue,  consisting  of  hydrogen, ■ oxygen,  nitrogen,  carbon, . sulpher 
and  phosphorous,  extending  from  x ■ o to  x ■ + ® where  x is  the  depth  of 
the  tissue.  It  is  assumed  to  be  placed  against  the  semi-infinite  non- 
capturing isotropically  scattering  reflector,  ■ which  extends  from  x * - » 
to  x * o.  The  source  consists  of  neutrons  entering  the  reflector  at 
x = - ® and  will  be  normalised  so  that  the  current  is  l/'JT  • It  is 
possible  to  show  that  in  this  case  the  flux  at  a detector  at  x * o in  the 
absence  of  the  tissue  is  unity.  The  reason  for  normalisation  to  unit 
flux  rather  than  to  unit  current  at  x ■ o is  that  flux  is  the  quantity 
measured  by  a detector.  Thus,  this  method  allows  of  comparison  of 
results  for  the  isotropic  source  with  those  previously  obtained  for 
neutrons  entering  the  tissue  at  right  angles  to, it.  The  numerical 
calculations  are  carried  out  in  this  report  for"  two  values  for  the  incident 
energy  namely  for  1 Mev  and  5 Mev. 

The  law  for  the  scattering  of  energy  flux  can  be  obtained  by  assuming 
that  the  nuclei  of  the  atoms  in  the  tissue,  other  than  the  hydrogen  atoms, 
are  infinitely  heavy  compared  with  the  mass  of  a neutron.  This  law  may  be 
described  fairly  accurately  as  a fourth  order  polynominal  in  the  cosine  of 
the  scattering  angle.  The  angular  distribution  of  energy  flux  can  then  be 
calculated  using  the  Spherical  Harmonics  Py  approximation.  This  process 
introduces  four  constants  of  integration  for  the  solution  in  the  tissue 
and  four  for  the  solution  in  the  reflector.  These  can  be  obtained  by 
equating  the  angular  distribution  in  the  tissue  to  that  in  the  reflector 
at  x = o 


2.  Equations  for  neutron  energy  flux. 

It  has  been  shown  that  the  biological  effect  at  a point  is  approximately 
proportional  to  the  energy  flux.  Let  this  be  denoted  by  e0(x)  where  x is 
the  depth  of  penetration  of  the  neutron  into  the  tissue  and  u is  the 
component  of  the  velocity  direction  vector  in  the  x direction. 

Let  us  now  consider  the  equation  for  the  neutron  energy  flux  in  the 
tissue.  This  is  given  in  T/R  273  where  it  has  been  shown  that  the 
angular  distribution  of  neutron  energy  flux  e(x,u)  satisfies  the  equation:- 


u 9e(x.  b)  + e(x,  u) 
Ox 


= _1 .//  e(x, u )[  JHi  + h(uo! 
2n  2 
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where  a neutron  travelling  in  direction^  is  scattered  to  a direction 

— , Pc  = J2  ‘ JL  , and  e0(x)  = /_  e(x,u)  du.  h is  the  probability  of 
collision  with  a proton,'  and  1— h the  probability  of  collision  with  a 
heavy  nucleus. 
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and  e(x,u) 


EjtxlPjfu)  . 


Substituting  these  expressions  into  1)  and  equating  coefficients  of  Pi(u) 
on  both  sides  of  the  equation  we  have 


e » ( go  ~ 1)  e0 


and  ( i +l)e,  et _ ^ 
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for  t > o. 


Since  we  are  using  a F?  approximation,  e is  neglected  yielding  the 
equations:-  8 
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The  solution  of  such  a set  of  equations  is  discussed  in  detail  in  T/R  231 
and  much  of  the  numerical  work  used  for  th^t  report  was  incorporated  in 
this  one. 

An  expression  is  obtained  for  enC x ) of  the  form 


* 

e„(x/  i)  = E an(vs)  e _V*X/  n,  * o,  1,  . . . .7. 

3=1 


The  ratios  an/a0  and  the  values  of  the  v'3s  were  obtained  from  previous 
calculations.  The  four  a0's  are  the  constants  of  integration  mentioned 
previously. 

The  equate  ns  for  the  energy  flux  in  the  reflector  have  now  to  be 
considered.  In  this  medium  the  angular  distribution  satisfies  the 
equation  for  isotropic  scattering  in  one  group  theory. 

9e(x,  p)  , , , 

i.e.  u + e(  x,  u>  * ? Eo(  x). 

0X 
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In  the  {^approximation  e(x,  u)  is  expanded  as  before 


Li  +1 

e(x,u)  * 2 — 6}<x)  PjCu)  . 

i-o  2 


The  4j(x)  in  this  case  satisfy  a set  of  equations  similar  to  2)  with 
go  " 1 and  gi  « o for  i >,  1 

Hence  it  can  be  shown  that 
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and  for  7 >,  k » 2 ek(x) 


The  formulae  for  the  Gk's  in  terms  of  v are  given  in  MDDC  236  and  may  be 
calculated  from  there. 

At  the  boundary  x » o,  e(x,  (i)  is  continuous.  Equating  coefficients 
of  the  Flip.)  , eight  equations  are  obtained  in  the  eight  unknown 
quantities  aoWJ,  ao(v2),  a0(vs),  ao(vJ  At  , Aa  , A*  , ■ A*  . 


3.  Numerical  Results. 

Two  values  of  incident  energy  were  considered,  namely  1 Mev  and  5 Mev. 
The  values  of  h Snd  of  the  mean  free  path  for  neutrons  in  soft  tissue  at 
these  energies  were  given  in  T/R  273  and  are  as  follows:- 


E in  Mev 

h 

liraean  free  path)  in  ems 

1 

•77  5 

2-82 

5 

•688 

6-37 

These  values  were  then  substituted  in  equations  2)  and  the  energy  flux 
eo(x/*)  determined  at  any  point  distance  x.  Now  the  energy  absorbed  per 
cc  of  tissue  per  second,  for  our  normalisation  is 

E e0(x/J) 
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To  convert  this  into  a dose,  measured  in  rep.  received  during  an  eight  hour 
period,  this  quantity  is  multiplied  by  a factor  10”4  x 4.9494. 

This  dose  is  plotted  in  Fig.  I,  as  a function  of  distance  into  the 
tissue,  for  the  cases  of  1 Mev  and  5 Mev.  Also  plotted  in  the  same 
diagram  are  the  corresponding  curves  for  neutrons  entering  at  right  angles 
to  the  surface. 

The  angular  distribution  of  emerging  neutron  energy  flux  for  the  case 
when  neutrons  come  from  the  reflector  is  plotted  in  Fig.  II. 

As  was  expected,  the  surface  dose  was  greater  for  the  isotropic  case 
than  that  when  the  neutrons  enter  at  right  angles  to  the  surface,  but  the 
dose  decreases  more  rapidly  for  greater  depths  in  the  tissue.  It  was 
found  that  the  isotropic  source  gave  increases  of  about  18°yfc  for  both  energy 
levels  on  the  surface  doses  for  the  albedo  case.  The  dose  in  each  case 
is  measured  as  the  ratio  of  the  actual  dose  to  the  flux  measured  at  a 
detector  in  absence  of  the  tisfeue. 

The  higher  surface  dose  is  due  to  the  fact  that  neutrons  entering  at 
directions  other  than  normal  to  the  surface  do  not  have  to  undergo  so  many 
collisions  in  order  to  be  reflected  back  to  the  surface.  Furthemore 
neutrons  leaving  the  tissue  can  be  reflected  back  into  it. 

It  should  be  noted  that  the  whole  of  the  incident  energy  is  absorbed 
in  the  present  case,  whereas  in  the  albedo  case  a small  fraction  is 
reflected  back.  It  will  be  remembered  that  the  current  is  l//~5  while 
in  the  albedo  case  it  is  unity  in  the  absence  of  the  tissue.  This  explains 
why  the  area  under  the  energy  deposition/distance  curve  for  the  isotropic 
case  is  about  l/vT  cf  that  under*  the  other. 

Tolerance  Flux  of  Fast  Neutron s . 

Assuming  a high  value,  go,  for  h the  relative  biological  efficiency 
factor  for  protons,  the  maximum  flux  of  fast  neutrons  from  a reflector, 
measured  at  a detector  in  absence  of  the  tissue,  which  in  8 hours  gives  a 
maximum  biological  effect  of  .lr  of  Y radiation  is  55  for  1 Mev  and  2« 
for  5 Mev. 

An  appreciable  proportion  of  the  surface  dose  is  due  to  neutrons  which 
have  been  reflected  to  and  fro  between  the  reflector  and  the  tissue.  If 
the  tissue  is  placed  in  direct  contact  with  the  free  surface  along  thc- 
whole  of  its  length  it  will  abserb  all  of  these  neutrons.  As  the  body 
is  moved  away  from  the  reflector,  however,  the  chance  of  the  tissue 
receiving  neutrons  which  have  been  reflected  in  this  way  is  diminished. 

The  past  of  the  surface  dcse  due  to  neutrons  which  are  coming  from 
the  reflector  after  undergoing  one  or  mere  reflc-cticns  between  it  and  the 
tissue  is 


/ e(o,u)  du  - 1 


This  is  about  10cyfc  of  the  whole  surface  dose. 

Therefore  if  the  body  is  at  a distance  from  the  reflector  which  is 
large  compared  with  the  dimensions  of  the  body  the  surface  dose  will  be 
about  lO'ft  smaller  than  indicated  in  Fig.  I. 
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FIG.  2 ANGULAR  DISTRIBUTION  OF  ENERGY  FLUX  AT 

THE  SURFACE  OF  THE  ISOTROPICALLY  SCATTERING 
MEDIUM  IN  THE  PRESENCE  OF  THE  TISSUE. 
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